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Abstract 

Motivated by gauge theory/string theory correspondence, a lower bound on the bulk 
viscosity of strongly coupled gauge theory plasma was proposed in [1]. We consider 
strongly coupled A/" = 4 supersymmetric Yang-Mills plasma compactified on a two- 
manifold of constant curvature /3. We show that the effective (1 -|- l)-dimensional 
hydrodynamic description of the system is governed by the bulk viscosity violating the 
bound of [1], once /3 < 0. 
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1 Introduction 

Relativistic hydrodynamics is an effective description of a nearly-equilibrium system at 
time (length) scales that are much longer (larger) than any characteristic microscopic 
scale of a system. In the simplest case (with no conserved charges), the dynamics of 
the hydrodynamic fluctuations in the system on a (p-|- 1) -dimensional manifold M.p+i 
with the metric 

ds^ ^Gf,^dx^'dx\ (1.1) 
is governed by conservation of the stress-energy tensor T^^, 

V^T^"' = 0. (1.2) 

The stress-energy tensor includes both and equilibrium part (with local energy density 
E and pressure P) and a dissipative part 11'''^, 

T^"" = SuV + PA'"' + W where A'''' = G^'' + uV . (1.3) 

Above, u'^ is the local (p-l-l)-velocity of the fluid with u^Uf^ — —1. Further, U.^'^Ui, — 0. 
In phenomenological hydrodynamics, the dissipative term 11'*'^, 

00 

n^^^^nr, (1.4) 

n=l 
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can be represented as an infinite series expansion in velocity gradients (and curvatures, 
for a fluid in a curved background), with the coeflicicnts of the expansion commonly 
referred to as transport coefficients. In (1.4) the subscript „ denotes the total number 
of the velocity da^ ■ ■ ■ daiUf} and/or the background metric da^ ■ ■ ■ da^Ga-y derivatives. 
The familiar example of the Navier-Stokes equations is obtained by truncating U^" at 
linear order in this expansion 

n'^'^ = nf (?], C) = -Tja^"" -CA^"'W-u, (1.5) 

where (for p > 1) 

^ 2V<'^M^> = A'^^A'^'^ (V«M/3 + V^Ma) - ^A'^^ (A^'^VaM/j) . (1.6) 

Notice that at this order in the hydrodynamic approximation we need to introduce only 
two transport coefficients, namely the shear^ rj and bulk ( viscosities. At the second 
order in the derivative expansion, Ilg^, one needs to introduce flve [2] or thirteen [3] 
additional transport coefficients (depending on whether or not the system is conformal). 

Holographic gauge theory/string theory correspondence [4] provides a useful general 
guidance about hydrodynamic transport coefficients. Thus, in [5], partly motivated by 
this correspondence, the authors (KSS) proposed a bound on the ratio of the shear 
viscosity to the entropy ratio^ 

s Att 

The KSS bound is obeyed in all known fluids in Nature. It is either saturated or satisfled 
in all explicit or phenomenological examples of holographic gauge theory plasma at 
inflnite coupling [5-9]. While the bound survives the leading order 't Hooft coupling 
corrections in four-dimensional holographic conformal models [10-12], it can be violated 
in strongly coupled conformal gauge theories with fundamental matter^ [13-15]^. 

In [1] we proposed a bound on the bulk viscosity^ of (inflnitely) strongly coupled 
gauge theory plasma 



>2 --C , (1.8) 

7] \p ' 



^To define the shear viscosity one needs p > 1. 
^We use h = ks = I- 

^More generally, in four-dimensional conformal theories with different central charges: a. 

'^See [16-20] for further discussion of the shear viscosity bound. 

^See [21-23] for the bound (and it violation) of a particular second-order transport coefficient 
"the effective relaxation time" . 
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or equivalently (using the universality of the holographic shear viscosity), 




(1.9) 



where is the speed of sound waves in plasma. In what follows we use the second 
version of the bound, eq. (1.9), as it is applicable even for p = 1. The bulk viscosity 
bound (1.9) is saturated in toroidal compactifications of conformal theories [1,24]. It 
is satisfied in various examples of string theory embedding of the holographic gauge 
theory/gravity correspondence: the strongly coupled J\f = 2* gauge theory plasma 
[25,26], the cascading gauge theory plasma [27,28], the mass-deformed Af — 4 SYM 
plasma with a non- vanishing U{1) R-charge chemical potential [29]. The bulk viscosity 
bound (1.9) is further satisfied in some phenomenological models of the holographic 
correspondence [30,31], however, it is violated in some other phenomenological model 



In this paper we address the question whether the violation of the bulk viscosity 
bound (1.9) is limited to the phenomenological models of the holographic gauge/gravity 
correspondence. A natural first guess is to replicate the strategy in the study of the KSS 
shear viscosity bound (1.7), i.e., to consider the effect of the gauge theory finite coupling 
and non-planar corrections on the purported bound. Unfortunately, such an approach 
can not lead to a reliable conclusion. Indeed, when the bulk viscosity originates from 
toroidal compactifications of conformal hydrodynamics (in which case the string theory 
dual corrections to finite coupling/non-planar corrections can be under control), the 
general arguments in [1] guarantee what the bulk viscosity bound continues to be 
saturated. On the other hand, a five-dimensional gravitational description of a generic 
non-conformal holographic correspondence involves scalar fields that originate from 
3-form RR fluxes of type IIB supergravity — in this case the full set of the higher 
derivative corrections to the supergravity is unknown. Clearly, we need to focus on 
explicit string theory examples of the holographic correspondence in the supergravity 
approximation. 

Our strategy in exploring the bound (1.9) is to consider compactifications of the 
strongly coupled conformal hydrodynamics on /c- dimensional spatial manifolds A4k, 
with /c < (p — 1), of constant curvature /3 7^ 0. Specifically, we consider hydrodynamic 
transport of strongly coupled N — A SU{N) SYM plasma on M.2, 



[32,33]. 




(1.10) 
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where is a round two-sphere, and is a hyperbohc space and G is a discrete 
subgroup of its SL{2, R) symmetry group. The quotient is assumed to be smooth 
and compact. L is a "radius" of Ai2- Since we are left with a single infinitely large 
spatial direction, the low-energy effective description of the compactification is given by 
(1 -|- l)-dimensional hydrodynamics. Notice that the latter hydrodynamic description 
can not be obtained from the compactification of the curved-space (3 + l)-dimensional 
hydrodynamics of = 4 plasma — the gradients of the background metric along 
the AI2 directions are generically^ large in the hydrodynamic limit, which invalidates 
the gradient expansion (1.4). Running ahead, we find that to leading order in ^ <^ 
1, the Al2-compactified hydrodynamics is affected by the third-order A/" = 4 SYM 
hydrodynamics on AI4 = i?^'^ x A^2- To first-order in the velocity gradients, the 
effective (1 -|- l)-dimensional hydrodynamic description of the theory is determined by 
a single transport coefficient — the bulk viscosity (. We extract ( from the dispersion 
relation of the sound waves in A/" = 4 plasma on R^'^ x AI2 propagating along the 
single uncompactified spatial direction, which we denote as '2;', 

to = cU-^^- q' + 0{q'), (1.11) 
s 

where 

and s is the entropy density. The speed of sound waves is determined from the equation 
of state as follows 



2 '-'-^ zz 



(1.13) 

/3=const 

where £ is the energy density, and Pzz = T^z is the equilibrium pressure in the plasma 
in the 2;-direction. 

The rest of the paper is organized as follows. In section 2 we discuss the regular 
black hole solution in AdSr, with the asymptotic boundary metric R^'^ x M.2, dual to 
an equihbrium thermal state of = 4 SYM plasma compactified on M.2- We compute 
the one-point correlation function of the boundary stress-energy tensor and find the 
general expression for the speed of the sound waves (1.13). For general values of 
the background geometry (and its thermodynamics) is found numerically. We present 
analytical results for the thermodynamics to leading order in ^ <^ 1. In section 3 we 



characteristic dimensionless parameter here is while a typical hydrodynamic parameter is 



f or M. 
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compute the dispersion relation of the sound channel quasinormal modes [34] in the 
black hole geometry of section 2, and extract the speed of the sound waves as well as 
the sound waves attenuation coefficient F = tt^. We present analytic results to leading 
order at high temperatures, 7^ <^ 1, and numerical results for generic values of In 
section 4 we discuss the bulk viscosity bound and its violation in = 4 SYM plasma 
on R^'^ X A42- We conclude in section 5. 



2 Thermodynamics of TV = 4 SYM plasma on i?^'^ x M2 

Effective five-dimensional gravitational action describing A/" = 4 SU (N) SYM on A^4 = 
R^'^ X A^2 in the planar limit {qym — )■ 0, — )■ 00 with A = Qym-^ — const), and for 
an infinite 't Hooft coupling A — )■ 00, is given by 

S^ = T7^I d'CV^iR + 12), (2.1) 

IDTTLrs J^g^ dMB=M4, 

where without loss of generality we normalized the radius of curvature of asymptotically 
AdS^ geometry to 1. Such a normalization implies 

= At . (2.2) 
5 27V2 V ; 

The background geometry dual to a state of the theory with translational invariance 
along the 2;-direction and the symmetry of is given by 

dsl = -cl dt^ + C2 I {dM2f + cl dz^ + cl dr^ , (2.3) 

where Cj = Cj(r) and {dM.2y is the standard metric on M.2- Given (2.3) we find the 
following second order equations of motion 

0=ci'-M£|)! + e;f^-^ + ^)-^-2c,c^ + ^%, (2.4) 
3c2 \3c2 C4 6C3 J 0C2C3 6 C2 



3C2 



0=4'-<f;HJ-ic;(^-^) + (l^-^)4-2c2c3 + f^#, (2.6) 

3C2 3 \C2Cx Cl/ \3C2 C4^J 6 C2 

as well as the first order constraint 

+ (c'2)^ + ^ - ^cl - &clcl + 2^ . (2.7) 
Cl C3C1 2 C3 

We explicitly verified that (2.7) is consistent with (2.4)-(2.6). 
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2.1 Asymptotics of the background geometry dual to a thermal state of 
the theory 

To describe a thermal state of strongly coupled J\f — A SYM plasma on R^'^ x M.2 we 
find it convenient to use the radial coordinate 

x = l--, (2.8) 

and further introduce^ 

, , a(x)q(x) , , a(x) 

Near the boundary, i.e., as a; — )■ 0+, the asymptotics of {a,g} are given by 



2 



Inx^ a; + O (a;-^/^) 



(2.10) 



They are characterized by two independent parameters: {/i , 02,0}- As we will see later, 
II is related to the temperature of the thermal state, and 02,0 determines the one-point 
correlation function of the boundary stress-energy tensor at thermal equilibrium. Near 
the regular horizon, i.e., as y = (1 — x) — > 0+, the asymptotics of {a , g} are given by 

a ^Ja^, + \4 iS2g'Mr - l) + 0{y')) , 

^ n ^ (2.11) 

A* 

They are characterized by three independent parameters: {oq > fl'o > }■ Notice that, 
given II, there are precisely four parameters 

{02,0 ,4,9o, 9i} , 

which is necessary to uniquely specify a solution for the second-order ODEs for {a, g}. 

Given (2.11) we can compute the Hawking temperature T and the Bekenstein- 
Hawking entropy density s of the black hole geometry (2.3) 

'2^^^'_ 1 ^ _ («o)'^o (2.12) 



II J 8(7oV 



^Thc two second-order equations of motion for {a,g} can be obtained from (2.4)-(2.7). As they 
are not very illuminated, we do not present them here. 
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For general the equations of motion for {a , g} are too complicated to be solved 
analytically — we solve them numerically. Analytical solution is possible though 
to leading order in the dimensionless parameter Using the high-temperature 
par ametrization , 

a{x) = (l + + ^ (^)) , 9{x) = ^ + j2 9ii^) + ^ (fj) ' ^2-^2) 

we find 

1 + 2x - N , 2x-x^- 2y/2x - 



diix) — — ; — - — %r flnfl — x) + arctanhVSx — + 
^ ^ 24{2x -x^)\ ^ ' ) 4&{2x - x^) 

\ 1 

gi{x) = — -arctanh\/2a; — x'^ — - ln(l — x) . 

Given (2.14), we find« 



(2.14) 



^'-^ 8^''%2+^(^^4y ' ^1-32+ 256 /i^+^U' 



(2.15) 



which implies (see (2.12)) 



2.2 Holographic renormalization and the boundary stress-energy tensor 

Holographic renormalization of strongly coupled A/" = 4 SYM plasma has been exten- 
sively discussed in the literature, see [35] for example. To render (dual gravitational) 
correlation functions finite, the effective action (2.1) has to be supplemented with the 
following set of counterterms: 

Set = I ^il + \P- ^AP'^Pki - P') In el , (2.17) 



SttGs JM,,e=c-^ V 2 12 



^3 



where 7 is the metric (2.3) restricted to C3 = e, and 



P = 7'''4- , Pij = \ (Pij - ^^7., ) ■ (2.18) 



^Wc get 0{(3/fi^) coefficient of g'l we actuaUy need to compute the second order correction, g2{x) 
in (2.13). Tliis is a straigfitforward extension and we omit tlie details. 
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In practice, we use expressions for tlie regularized one-point correlation function of 
the stress-energy tensor obtained in [36], in particular, see eq. (3.52), in the conformal 
limit. We find: 

8.G.r„=^£=2/-2/a,. + l/.M„^, ^^^^^ 

8XG5 r„ P„ = i / + 2/ 02,„ - 1 /S^ In ^ , 

where A is an arbitrary (fixed) renormalization scale associated with the ambiguity^ of 
defining the stress-energy tensor of a theory on curved background manifold A44 with 
broken supersymmetry. 

Given (2.19), it is straightforward to compute the speed of the sound waves (1.13) 
in Al2-compactified hydrodynamics oi J\f — A SYM plasma: 



d£ 



_ 1 + 4a2,o(/3) - 2/3 „(/3) - p_ 
;3=const 3 - 4a2,o(/3) + 2/3 a',M + ^Z?' ' ' 



where the prime denotes derivative with respect to /3. In section 3 we compare the 
thermodynamic prediction (2.20) for the speed of the sound waves with direct compu- 
tations from the dispersion relation of the sound channel quasinormal modes. 
Using the asymptotic expansions (2.15) and (2.16) we find 



(2.21) 



andi° 

2 



3 18//2 432/^4 ^ \^^6 



ci=o-T^.-J^.+0 . (2.22) 



Notice from (2.16) that since 



2 V 4772^2 

the basic thermodynamic relations 

-P,, = T = £-Ts, dS^Tds, (2.24) 

are (analytically) satisfied to order O (^j^^ ■ 

^The same ambiguity appears in computation of the thermal stress-energy tensor of strongly cou- 
pled A/" = 2* plasma [37]. 

^°To compute to order 0{/3^/iJ^) we do not need the 0{j3'^ / ji^) coefficient in the expansion of 
a2,o, see (2.15). 
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3 Sound channel quasinormal mode of AdS^ black hole with 
i?^'^ X At 2 asymptotic boundary 



Following [34], the dispersion relation for the sound waves in A^2-compactified A/" = 4 
SYM plasma is identified with the dispersion relation for the sound channel quasinormal 
modes propagating in the 2;-direction of the black hole geometry (2.3). We briefly 
outline the construction of the corresponding fluctuations. 
Consider the following decoupled set of metric fluctuations 

hu(t,z,r) ^cI{t)Hu{t) e--*+^«^^ 

ht,{t, z, r) =cl{r) HUr) e--*+^«^^ , (3.1) 
h,,{tz,r) =cl{r) Hss{r) 6"'-*+^^^^ 
hUt,z,r) =cl{r)HUr) e— *+^<'^^ 

where 

^if'^dede = ^{dM2r. (3.2) 

From the equations of motion for the fluctuations (3.1) we obtain 4 second-order linear 
ODEs for {Hft, Hfz, H^z, Hgs} as well as 3 linear flrst-order constraints associated with 
the diffeomorphism-flxing conditions 

h^r — hzr ~ hff — . (3.3) 

The combination 

^ + 2 - 2 f ^ - 4^ V 2^ ^?^3 Hu , (3.4) 

tt) \C2C2 tt)^c2C2/ 

is invariant under the residual diffeomorphisms (for the gauge fixing (3.3)) and satisfies 
the following decoupled linear ODE^^ 



{ci,2,3(r)},m,q 



(3.5) 



— Z'^ + Ci Z'jj + C2 Zh , Ci — C, 
In terms of the radial coordinate (2.8), the hydrodynamic limit takes form 
Zh{x) = (l-x)-^" {^ZHfl{x)+n ZH,i{x)+0{q^)^ , m = c, q-iP q'+e'(q=^) , (3.6) 

^^Thc expressions for Ci are too long to be presented here — they are available from the author 
upon request. 
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Figure 1: (Colour online) The speed of the sound waves in A^2-compactified A/" = 4 
SYM plasma as a function of j2§T)^- "^^^ dashed blue curves indicate the speed of the 
sound obtained from the quasinormal mode dispersion relation, see (3.6). The solid red 
curves indicate the prediction for the speed of the sound waves from the equilibrium 
thermodynamics, see (2.20). 



with the following boundary conditions 

lim Zh^q = 1 , hm Zh,i = , 

ZH,i = 0{x) , as X — > 0+ . 

We omit the details of solving the boundary value problem (3.5)-(3.7) for {zh,o, zh,i} 
and present only the results^^. First of all, we find 



± e. ^ 1 _ I 4 _ 0,004009,4) ^ + O (^ 

V3 36 /i^ Va* 

r=-H — + oi — 

3 18 /i2^ V/"^ 



(3.^ 



Notice that the speed of the sound waves is in perfect agreement with the predictions 
from the thermodynamics (2.22) — the coefficient of ^ agrees with an accuracy of 



3j 

oc 10-". 



Figure 1 presents the speed of the sound waves in A^2-compactified M = 4 SYM 
plasma as a function of j2§r)^- '^^^ dashed blue curves are obtained from the dispersion 
relation of the sound quasinormal modes (3.6), while the solid red curves indicate the 
prediction for the speed of the sound waves from the equilibrium thermodynamics, see 



^For more details on the solution procedure see [26]. 
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Figure 2: (Colour online) The attenuation of the sound waves (solid blue lines), F = vr^, 
in A^2-compactified TV = 4 SYM plasma as a function of j^fy- The solid green line 
indicates the high-temperature prediction (3.8). The vertical dashed red line indicates 
the critical temperature at which the speed of the sound waves vanishes. 



(2.20). Notice that at 



^ - 0.570580(8) , cl 



= 0, (3.9) 

r=T„ 



the speed of the sound waves squared vanishes, and continuous to the unstable branch 
with < 0. In the vicinity of the critical point oc (l — ^'Y^'^ — same critical 
phenomena has been observed in = 2* gauge theory plasma [37], and the cascading 
gauge theory plasma [28]. 

Figure 2 presents the attenuation F = vr^ of the sound waves in A^2-compactified 
A/" = 4 SYM plasma as a function of ^^j^^a (solid blue curves). The solid green line 
indicates the high-temperature prediction (3.8). The vertical dashed red line indicates 
the critical temperature Tc, see (3.9). Notice that as in the case of M = 2* gauge 
theory plasma [1], the bulk viscosity remains finite at the critical point. 



4 Bulk viscosity bound and its violation in A/^ = 4 SYM plasma 

on i?^'^ X M2 

The bulk viscosity bound proposed in [1] , in the case of strongly coupled A/" = 4 SYM 
plasma compactified on A^2 (see (1.10)) reads 

- = ^{^-cl), A>1. (4.1) 
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Figure 3: (Colour online) Violation of the bulk viscosity bound A > 1 (see (4.1)) in 
= 4 SYM plasma compactified on 2 as a function of the conformal symmetry 



Clearly, if /3 < 0, i.e., strongly coupled A/" = 4 SYM plasma is compactified on S^, 
the bulk viscosity bound is violated. The leading (at high-temperature or equivalently 
small A^2 curvature compactification) violation of the bulk viscosity bound is oc /3 oc -R. 
It comes from 2-derivatives of the metric along A^2 directions. When viewed from the 
perspective of higher-order hydrodynamics of A/" = 4 SYM plasma on A^4 = i?^'^ x A^2, 
such a violation comes from the third-order dissipative term, Ilg'^. 

Figure 3 presents parameter A in the bulk viscosity bound (4.1) as a function 
of the conformal symmetry breaking parameter (| — c^) of strongly coupled J\f = A 
SYM plasma compactified on A^2- Notice that A < 1 always, as long as /3 < 0, 
i.e., the compactification manifold A^2 is a constant curvature higher genus Riemann 
surface. The violation is rather strong at low temperatures: it is ~ 51% at the lowest 
temperature we accessed numerically 



breaking parameter ( 



1 
3 




From (3.8) we find 




(4.2) 



mm 



(27rT)2 



-30.379. 



(4.3) 
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5 Conclusion 



In this paper we constructed a specific string-theoretic counter-example to the bulk vis- 
cosity bound in (infinitely) strongly coupled gauge theory plasma proposed in [1]. We 
observed that compactification of the higher-dimensional hydrodynamics on curved 
manifolds results in Navier-Stokes (first-order) hydrodynamics with transport coef- 
ficients that are sensitive to higher-order dissipative terms of the higher-dimensional 
hydrodynamics. In the small-curvature limit of compactifications, the sensitivity starts 
with the third-order dissipative terms of the higher-dimensional conformal hydrody- 
namics. Since fiat-space compactifications of the conformal hydrodynamics saturate 
the bulk viscosity bound of the effective lower- dimensional hydrodynamics [1], confor- 
mal hydrodynamics compactifications on curved manifolds are guaranteed to violate 
the bound for the judicious choice of the compactification manifold curvature. In a 
specific example of strongly coupled A/" = 4 SYM plasma we showed that violation 
occurs for compactifications on negative curvature two-manifolds. 

Since A/" = 4 SYM contains conformally coupled scalars, one wonders whether the 
violation of the bulk viscosity bound can be attributed to the presence of tachyons"^^ 
in the theory, when compactified on S^. We do not believe this to be the case: in 
J\f = 2* gauge theory plasma [37] one can study mass deformations with ^ < 0; such 
deformations do not violate the bound [1] since, for example at high-temperatures, the 
bulk viscosity is affected at order j^. 

We did not explore in details the issue of stability of A^2-compactifications of A/" = 4 
SYM discussed here. We did verify though that, at least at high temperature, there 
are no instabilities of metric fluctuations whose wave-functions are constant over A^2- 
Likewise, there are no instabilities of minimally coupled scalar (for example, a dilaton) 
in the background geometry (2.3), again, provided its wave-function is constant over 
A^2- In case of compactiflcations (which lead to the violation of the bulk viscosity 
bound) fluctuations with non-trivial wave-function on might be projected by the 
action of G on the quotient = H^/G. The latter fact prevents making a general 
statement about the stability of the compactifications. 

Finally, in this paper we considered compactifications of A/" = 4 SYM on A/l2 
which completely break the supersymmetry. It would be interesting to extend analysis 
presented here to supersymmetric (twisted) compactifications of A/" = 4 SYM discussed 

^^Tachyons do not immediately imply the instabilities in the theory. 
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in [38]. 
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